Future CMB ISW-Lensing bispectrum constraints on modified gravity in the 
Parameterized Post-Friedmann formalism 



(N 

o 

(N 

> 
O 

(N 



Bin Hu 1 , Michele Liguori 1,2 , Nicola Bartolo 1,2 , and Sabino Matarrese 1,2 
1 INFN, Sezione di Padova, Via Marzolo 8, 35131 Padova, Italy 
2 Dipartimento di Fisica e Astronomia "G. Galilei", 
Universitd degli Studi di Padova , Via Marzolo 8, 35131 Padova, Italy 
(Dated: November 22, 2012) 

We forecast the constraints on both Hu-Sawicki model and Bertschinger-Zukin model of mod- 
ified gravity within the Parameterized Post-Friedmann (PPF) formalism for the Planck satellite 
experiment by performing the joint analysis of ISW-Lensing bispectrum and CMB power spectrum. 
We find that, even considering the temperature-temperature mode of CMB power spectrum only, 
Planck data are expected to reduce the error bars on the modified gravity parameter Bo (related 
to the present value of Compton wavelength of the extra scalar degree of freedom) at least one 
order magnitude compared with WMAP. The spectrum-bispectrum joint analysis can further im- 
prove the results by a factor ranging from 1.14 to 5.32 depending on the specific modified gravity 
model. One of our main results is that the cross-correlation between ISW-Lensing bispectrum and 
power spectrum can be safely neglected when performing the joint analysis. For simplicity, we only 
investigate the likelihood of one parameter (-Bo) and fix all other cosmological parameters to their 
best-fit values in WMAP7yr results. 
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I. INTRODUCTION 

Cosmic acceleration can arise from either an exotic 
form of energy with negative pressure, namely "dark en- 
ergy" , or a modification of gravity that appears on large 
scales. As shown in [1], every expansion history that can 
be parameterized by a dark energy model with pDE(lna) 
can be reproduced by a one parameter family of f(R) 
models, which approaches the Einstein-Hilbcrt action in 
the high curvature limit. Hence, the observations of the 
cosmic background expansion cannot help us to distin- 
guish the viable modified gravity from dark energy mod- 
els. However, the dynamical evolution of fluctuations can 
provide us a tool to test gravity theory (e.g. [2-14] and 
references therein). The main difference between mod- 
ified gravity and dark energy models is the anisotropic 
shear on large scales. In most of the dark energy mod- 
els the large scale anisotropic shear vanishes, while, in 
modified gravity models usually it does not. The deeper 
physical reason for this phenomenon is that in modified 
gravity usually the evolution of the Newtonian potential 
\1/ deviates from that of the spatial curvature potential 

A viable modified gravity model should at least satisfy 
the following three conditions: first, it should pass the 
stringent tests in the local solar-system; second, it should 
respect the standard expansion history and accelerate the 
late-time expansion; third, it should avoid modifying the 
physics at recombination, which has been investigated by 
various CMB experiments with a high precision. 

The studies of modified gravity models, in principle, 
can be classified into two different frameworks (see, e.g., 
the discussion in Rcf. [19]). The first is a model depen- 
dent method. One can start from a specific Lagrangian, 
investigating its dynamical behavior to finally give its 
predictions. Various viable modified gravity models have 
been proposed (see the review article [15]). In this pa- 



per we mainly focus on f(R) models (see [18] for a re- 
view), such as Starobinsky model [16], Hu-Sawicki (HS) 
model [17]. The other method is inspired by the Param- 
eterized Post-Newtonian (PPN) approach for the solar- 
system test. In this case one can try to build a more or 
less model independent framework, in which many mod- 
ified gravity models can be parameterized in an unified 
way. The simplest idea is to directly generalize Edding- 
ton parameter (7 = [20] in static solar-system 

to a one unknown function of space and time "f(t, x) 
in Friedmann universe. Many studies (e.g., [19, 21- 
25]) based on this framework show that this parame- 
terization works quite well with the large scale struc- 
ture data. However, as stated in [26] using a parame- 
terization with insufficient freedom significantly tightens 
the apparent theoretical constraints. Still in the same 
spirit of the standard solar-system PPN approach, an- 
other viable parameterization has been proposed, the so 
called Parameterized Post-Friedmann (PPF) framework 
[27-29], which allows to describe (and match) three dif- 
ferent regimes of gravity modifications to GR, namely 
the super-horizon regime, the (sub-horizon) quasi-static 
regime (with an intermediate regime between the two) 
and a non-linear regime on smaller scales. The current 
constraints on general f(R) models within the PPF for- 
malism are B < 0.42(95%C.L.) by using CMB and ISW- 
galaxy correlation data, and B n < 1.1 x 10~ 3 at 95% 
C.L. [31]. using a larger set of data, such as WMAP5, 
ACBAR, CBI, VSA, Union, SHOES, and BAO data. Be- 
sides what mentioned above, there exist many other pa- 
rameterizations [32-35] . 

On the other hand, from the theoretical point of view, 
it has become clear that the cross-correlation between the 
Integrated Sachs- Wolfe (ISW) effect [36] and the gravi- 
tational Weak Lensing (WL) in CMB anisotropics can 
provide a wealth of cosmological information, including 
the dynamics of cosmic acceleration [37-41]. Recently it 
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has been shown that this effect will provide the most im- 
portant contamination to the local-type primordial bis- 
pectrum [42-48]. In order to detect the primordial non- 
Gaussian signals, one needs to accurately subtract this 
secondary anisotropy in the Planck data. In fact vari- 
ous forecasts indicate that the Planck mission can mea- 
sure this secondary effect with a high precision [42-48]. 
From the observational point of view, determining the 
cosmological parameters will require not only additional 
measurements, but also characterizations of the data be- 
yond the power spectrum. ISW-Lensing bispectrum can 
represent one of these novel statistical variables. In this 
paper we focus on constraining modified gravity models 
via future Planck ISW-Lensing bispectrum data. 



II. PARAMETERIZATION OF MODIFIED 
GRAVITY 



The mathematical foundation of the validity of PPN 
formalism is Birkhoff's theorem, which provides an ex- 
tremely simple form of Schwarzschild solution in any geo- 
metric theory. However, compared with the PPN formal- 
ism, PPF is a parameterization of a scries of differential 
equations, rather than a specific static/stationary solu- 
tion. This is because what we are interested in cosmology 
is not only inhomogeneous but also time- varying quanti- 
ties. Fortunately, in some dominant regime for measures 
of the rate of structure growth, such as the weak lens- 
ing and peculiar velocity surveys, the time derivatives of 
perturbations can be neglected relative to spatial deriva- 
tives. This approximation leads to the so called quasi- 
static approximation, which can reduce the time differ- 
ential equations into simple algebra relations. However, 
on the super-horizon scales the long wavelength perturba- 
tions have enough time to feel the changes of the gravita- 
tional potentials, consequently, the time derivative terms 
in the differential equations cannot be neglected. One 
explicit example is the ISW effect. In this section, wc 
first recall one of the conventional PPF parametrization 
which is mainly inspired by the quasi-static dynamics, 
the Bertschinger-Zukin parameterization [19]; and then 
we introduce a more suitable parameterization for our 
purposes, the Hu-Sawicki PPF formalism, which prop- 
erly takes the super-horizon perturbation evolution into 
account, even for large value of Bq (/ho)- 



A. Bertschinger-Zukin parameterization 



governed by two modified Poisson-like equations 

-k 2 § = 4irGi(t,k)a 2 p m S , 
-k 2 ^> = 47rG*(i,k)a 2 p m (5 , 



(1) 
(2) 



with G*(i,k) = -f(t, k)G$(i, k) 1 . Given the above ob- 
servation, Bcrtschinger and Zukin proposed the following 
parameterization for f(R) model 



G, 



1 + aik 2 a s 



Go l + a 2 k 2 a s 



l + Pik 2 a s 

7 = IT^V' (3) 



with 



*(t,k)=7(t,k)$(t,k) 



(4) 



where Go is the Newton constant at the present epoch. 
We can easily see that in the quasi-static approxima- 
tion time derivative terms are neglected so that long- 
wavelength perturbations cannot be properly taken into 
account (e.g., 7 — > 1 for k — > 0, while in general, mod- 
ified gravity theories can show distinctive features also 
on the longwavelength regime). However, authors of [61] 
recently demonstrated that for small values of B (fm), 
the parameterization of [19] is practically good enough 
for the current data analysis purpose. In our following 
studies, we also investigate this kind of parameteriza- 
tion. Constraining modified gravity models via the ISW- 
Lensing bispectrum forces to take the ISW effect into 
account properly. Given this request, in what follows we 
turn to Hu-Sawicki PPF formalism in that it allows to 
join the superhorizon to the subhorizon regime. 



B. Hu-Sawicki model in PPF formula 

The PPF parameterization is defined by 3 functions 
<7(lna, kff), ./^(lna), /c(lna) and 1 parameter cr- They 
correspond to the metric ratio, the super-horizon rela- 
tionship between the metric and density, the deviation 
of Newton constant on super-horizon scale from that 
on quasi-static regime, and the relationship between the 
transition scale and the Hubble scale [27]. The 4th or- 
der nature of f(R) theories provides enough freedom to 
reproduce any cosmological background history by an ap- 
propriate choice of the f(R) function [1]. For simplicity, 
in this paper we concentrate on one specific parameter- 
ization form, namely Hu-Sawicki model [17], which can 
satisfy the local solar system tests 



f(R) 



-m 



ci(i?/m 2 ) n 
c 2 (R/m 2 ) n + l 



(5) 



As a straightforward extension of the PPN formal- 
ism, Bertschinger-Zukin formalism (BZ) [19] proposed 
a quite efficient parameterization of modified gravity in 
the quasi-static regime, where the dynamical behavior 
of the temporal and spatial gravitational potentials are 



Notice that the metric convention of Bcrtschingcr-Zukin is differ- 
ent w.r.t. Hu-Sawicki, we have: <I>bz = ^HS an d ^BZ = — ^HS- 
In this paper, except in this subsection, we follow Hu-Sawicki 
conventions. 
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where the mass scale reads 



m 2 ee H%fl m = (8315Mpc)" 



0.13 



(6) 



The non-linear terms in f(R) introduce higher-order 
derivatives into this theory. However, we are more famil- 
iar with a second-order derivative theory. Fortunately, 
we can reduce the derivatives to second order by defining 
an extra scalar field \ = df/dR, namely the "scalaron", 
which absorbs the higher derivatives. The Compton 
wavelength of the scalaron is defined as 



B = 



Irr R , H 
1 + fn H> ' 



(7) 



with f R = df/dR, f RR = d 2 f/dR 2 and ' ee d/dlna. 

At the high curvature regime, (6) can be expanded 
w.r.t. m 2 jR as 



lim f(R) 



Ci 2 

— m 

C2 



a 2 
3 m 



m 



(8) 



From Eq. (8) we can see that, the first and second terms 
represent a cosmological term and a deviation from it, 
respectively. In order to mimic a ACDM evolution on 
the background, the vale of Ci/o^ can be fixed as [17] 



(-2 



6 



(9) 



By using the relation (9) the number of free parame- 
ters can be reduced to two. From the above analysis, 
we can see that, strictly speaking, due to the appear- 
ances of correction terms to the cosmological constant, 
Hu-Sawicki model cannot exactly mimic ACDM. Since 
m 2 1 R increase very fast with time, the largest value (at 
the present epoch) is m 2 /R ~ 0.03, the largest deviation 
to ACDM background happens in the case n = 1 with 
1% errors, corresponding to m 2 /Rc2 ~ 0.01 in Eq.(8). 
For large n values, such as n = 4, 6 we can safely neglect 
this kind of theoretical errors. We have checked that for 
n = 1 case, this 1% deviation from ACDM will bring us 
a 10% errors in the variance of the parameter Bq, while 
for n — 4, 6 our results are not affected. 

Without loss of generality we can choose the two free 
parameters to be (n, c%). However, as shown in [1], for 
more general f(R) models the ACDM evolution on the 
background can be reproduced exactly by only introduc- 
ing one free parameter. This means that there exists 
some degeneracy between the two parameters, as shown 
in Fig. 1. Usually General Relativity (GR) is recovered 
when 



Bn = 



(10) 



Then, from Fig. 1 we can see no matter what value 
n takes, we are always allowed to set B Q = by ad- 
justing C2- Furthermore, in order to mimic ACDM on 
the background, C2 and n need to satisfy one constraint: 



the first term in the denominator of (5) should be much 
larger than the second. The blank area in Fig. 1 cor- 
responds those parameter regions that are ruled out by 
C2(R/m 2 ) n > 1. Actually, this condition gives 



nmax 
B 



0.1, (n = l), 
1.2, (n = 4), 
4.0 , (n = 6) . 



(11a) 



In the approximation (8), we assume m 2 / R —¥ 0, how- 
ever, with the expansion of the universe this approxima- 
tion starts to be not so accurate for our purposes. In 
order to obtain our results as accurate as possible, in the 
following calculation we use the exact expression (5) in 
place of the approximated one. 




FIG. 1. The degeneracy between parameters C2 and n. 

We have decided to study the Hu-Sawicki model within 
the PPF formalism for the particular relevance of this 
model, and to exploit the generality of the PPF for- 
malism for the study of the perturbations, using for the 
computations the publicy available PPF module [29] of 
C AMB [30] . In Ref. [49] we apply the same formalism to 
more general f(R) models and other models of modified 
gravity. As stated previously, the idea of the Hu-Sawicki 
PPF formalism is to parameterize gravity theory in order 
to join different scale regimes. In what follows, we will 
briefly review this parameterization. Firstly, we define 
the metric ratio as 
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$ + * 



(12) 



In principle, g is completely equivalent to the function 
j(t,k) introduced by Bertschinger and Zukin in [19]. 



1. Super- Horizon Regime 

On super-horizon scales, the dynamical equation of 
motion for the spatial curvature potential has the fol- 
lowing form [27] 

H" B' 



1 - 



ff_ 

IT 



H' 
H" 



1 — B 
B' 
1 — B 



4 1 * 



$ = 0, (fcff^O). (13) 
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Furthermore, the comoving curvature conservation rela- 
tion (C = 0) gives 



<f = 



l-B 



{k H -> 0). 



(14) 



On super-horizon scales, because the spatial gradients of 
the perturbations can be neglected compared with their 
temporal evolution, the metric ratio reduces to a function 
5sh (In a) dependent only on time 



g(]na,k H = 0) = gsH(lna) = 



$ + * 
$ - f 



(15) 



The relevant phenomenon about the decay/growth of the 
gravitational potentials in the CMB photons is the ISW 
effect. In Fig. 2 we compute the ISW effect in the mod- 
ified gravity models we consider in this paper. Our cal- 
culations demonstrate that ISW effect is quite sensitive 
to the value of the parameter n in the Hu-Sawicki model. 
The left plot shows that, for the Hu-Sawicki model with 
n = 1, the ISW effect is always decreasing with increas- 
ing of Bq. However, for n = 4 case, the ISW effect firstly 
decreases monotonically then bounce again after reach- 
ing some critical values. This behaviors are also reflected 
on our likelihood analysis (see later, Fig. 5). 




FIG. 2. ISW effect in modified gravity models: Bertschinger- 
Zukin (BZ) model, Eq.(3) and f(R) Hu-Sawicki model, 
Eq.(5). 



below the Compton wavelength. 

3. Matching the superhomzon and the subhorizon regimes 

In order to match the super-horizon scale behavior, an 
additional term V must be introduced to the modified 
Poisson equation [27] 



= AnGa p m A. r , 



$_ + r 

where T satisfy the following equation 
(1 + 4k%) \V + T + 4k%{T - / G $_ 



with the source term 

s 



3 H„ 



1 H 

-2g 



2H 2 w 



-(i + /c) 



s 



Vrn 

k H 



(18) 



(19) 



1 



(20) 



In Eq.(19), the coefficient cr represents the relationship 
between the transition scale and the Hubble scale, and 
the function the relationship between the metric and 
the density perturbation. For f(R) models, we have 



c r = 1 , f c =c c g, c c » 1/3 . 



(21) 



Up to now, the only unknown function in the intermeadi- 
ate regime is the metric ratio. An interpolating function 
is proposed in [27] as 



with 



/(In a, k) 



-1/3 , n g 



ffSH + ffQs(c g fcg)" g 

i + (c g k H y* 



(22) 



c g = 0.71 y/B{$ . (23) 



III. 



ISW-LENSING BISPECTRUM 



2. Quasi-Static Regime 

In the quasi-static regime, all the time derivative terms 
are neglected w.r.t. the spatial gradients. At the linear 
perturbation order, this brings a lot of simplifications, 
because we can reduce the ordinary differential equations 
into algebraic relations. As an example, the modified 
Poisson equation takes the following form [27] 

= J^ a 2 pmAm t ( 16 ) 

1 + JG 

where fa represents the deviation from the standard 
gravitational constant (here $_ = ($ — 4 r )/2). Further- 
more, for f(R) models, we have 

fa = fR, <?^<?qs = -1/3, (17) 



The primordial local non-Gaussianity is generated 
when the short-wavelength fluctuations are modulated 
by the long-wavelength modes which cross the horizon 
much earlier than the former. This kind of interactions 
will induce some correlations between the random fluctu- 
ations on large scales and those on small scales. In princi- 
ple, any mechanisms through which small and large scale 
fluctuations couple together will generate similar features 
as those in the primordial local-type non-Gaussianity. 

The late-time cross-correlation between ISW and Weak 
Lensing provides one such mechanism. On one hand, 
when the universe starts to accelerate at late-time, the 
decay of the gravitational potential generates the sec- 
ondary anisotropics in the CMB anisotropics, the ISW 
effect 

^(n)|isw - f d X ($ - *),r(A,x) , (24) 
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where r and x are the conformal time and comoving dis- 
tance, respectively. On the other hand, the gravitational 
lensing induced by the fluctuations of matter density de- 
flects the paths of CMB photons as they travel from the 
last-scattering surface to the observer 



ST(h) = <5T(n + , 

~<5T(n) + (d^)-(dST) (n) , 



where the lensing potential is defined as 



(n) = - 



X* 



X 



x*x 



(n, X 



(26) 



From Eq.(24) and (25), we can see that Weyl poten- 
tial ($ — Vf) sources both ISW and WL, so the long- 
wavelength mode from ISW couples with the short- 
wavelength mode from WL. Consequently, this kind 
of correlation will induce the secondary ISW-Lensing 
bispectrum, which has been recognized as the most 
dominant contamination of local-type primordial bispec- 
trum [42-47] 

Going to harmonic space and computing the reduced 
bispectrum, one can obtain [37, 38, 43, 59] 



^ISW-L 



-h(h+l) + h(l2 + l)+l 3 (k + l) 



+5perm. 



(27) 



jsw 



). Since the gravitational lensing 



where Cf T = ( 
is just remapping the primordial CMB temperature and 
polarization fields via the gradient of a foreground po- 
tential, from the above expression we can see that ISW- 
Lensing bispectrum includes information from both the 
primordial fluctuations and the secondary ones. This 
fact is important for the understanding of the spectrum- 
bispectrum joint analysis, because the bispectrum can 
correlate with the power spectrum through both the 
temperature-temperature mode of power spectrum (C ; T ) 
and the secondary lensing-temperature cross power spec- 
trum (Cf T ), as detailed in Sec. IV C. 

In Fig. 3, we plot the lensing-temperature cross power 
spectrum for several different gravity models. From this 
figure, we can explicitly see the non-linear dependence 
of Cf T on B : firstly, in the GR limit Bq -> (red 
solid curve), the cross-correlation is large and positive; 
secondly, with increasing of Bq, the cross-correlation de- 
creases and oscillates around zero (green dotted curve 
for n = 1, pink dashed curve for n = 4 Hu-Sawicki 
models, and black dotted-dashed curve for Bertschinger- 
Zukin model); finally, with further increasing of Bq, the 
cross-correlation amplitude increase again but negatively 
(dotted deep blue curve for n = 1, dashed-dotted light 
blue curve for n = 4 Hu-Sawicki models, and dashed- 
dotted orange curve for the Bertschinger-Zukin model). 
The reasons of these phenomena lie in the non-linear de- 
pendence of ISW effect on Bq, i.e. the first few multi- 
poles firstly decrease and then bounce with the increas- 
ing of Bq, see Fig. 2. The other important feature is 
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FIG. 3. The cross spectrum between lensing potential and 
temperature i 3 C* T . 




FIG. 4. The squeezed slice of reduced bispectrum I b4,i,i+4- 



that the cross-correlations are suppressed greatly after 
the first hundred multipoles, it almost vanishes in the 
regime I > 1000, because the cross-correlation is driven 
by the late time ISW effect. 

We also investigate the squeezed profile of reduced 
ISW-Lensing bispectrum, see Fig. 4. From this figure, we 
can easily see the modulation of short-wavelength mode 
by the long-wavelength mode. 



IV. FORECASTS AND RESULTS 

In this section, we will present the joint likelihood anal- 
ysis of the CMB power spectrum and ISW-Lensing bis- 
pectrum for several modified gravity models by using the 
specifications of WMAP and Planck experiments. In this 
paper, we implement the Hu-Sawicki model by using the 
PPF module [29] of CAMB [30] and the Bertschinger- 
Zukin model by using MGC AMB [24] . Our results show 
that the constraint on Bq is quite model dependent. 
However, generally speaking, the ISW-Lensing bispec- 
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TABLE I. Bo error bars (68% CL) from power spectrum \ 2 
for the Hu-Sawicki model 



n 


WMAP 


PLANCK 


1 


4.20 x 10~ 2 


6.68 x 10~ 5 


4 


6.14 x 10" 1 


2.06 x 10~ 3 


6 


4.79 x 10 _1 


1.51 x 10~ 2 



TABLE II. Bo error bars (68% CL) from power spectrum \ 2 
for the Bertschinger-Zukin model 



Ideal (Z max = 1000) 


PLANCK 


1.63 x 10 _1 


1.10 x 10~ 2 


a In order to compare with [41] 


we forecast Bo error bars in the 


ideal experiment with l max = 


1000 instead of WMAP. 



trum can constrain cosmological parameters at the same 
level as the power spectrum does. Consequently, the joint 
analysis of the power spectrum and the ISW-lensing bis- 
pectrum can improve the B error bars by a factor 1.14 
to 5.32. 



A. Power spectrum 

In this subsection, we forecast constraints on Bq from 
observations of the CMB temperature power spectrum. 
To this aim we firstly consider the simple likelihood ex- 
tracted from the chi-squared variable: 





i n=4, PLANCK 



0.005 0.01 0.015 



3000 



1=1 



c, 



fid 



(28) 



where C; is the lensed "TT" power spectrum, and the 
variance is 



4„ 



(21 + 1) 



(29) 



Ni and hi denoting the experimental noise and beam re- 
spectively (listed in Tab. VII). We allow only Bq to vary 
while fixing all other cosmological parameters to their 
WMAP7yr best- fit values, as in [41]. The B -likelihood 
obtained in this way is shown in figure Fig. 5 and 6. In 
agreement with previous studies [31, 41], we find that 
its behaviour is highly non-Gaussian and for this reason 
we cannot rely on a standard Fisher matrix approach in 
order to forecast error bars on B - Despite this issue 
we still computed a full 7-parameter Fisher matrix as a 
mean to get an order of magnitude estimate of the cor- 
relation between B and standard ACDM parameters. 
Our results are summarized in detail in Appendix A and 
show that there are no significant degeneracies between 
Bo and other parameters (the largest cross-correlation 
coefficient we find is n b h 2 -B Q = -0.18). This in turn 
justifies the simple 1-parameter likelihood approach used 
in this section in order to study Bq error bars. 

We investigate the likelihood function for both WMAP 
and Planck experiments in three Hu-Sawicki models with 
different values of n = 1,4,6 in Fig. 5 and 6. Further- 
more, in order to compare with the results in [41] we 
also include the Bersthinger-Zukin model in Fig. 6. The 
corresponding 68% CL's are summarized in Tab. I and 
II. From the above results, we firstly see that Planck is 



FIG. 5. WMAP/PLANCK experiments power spectrum like- 
lihood of Bo, for n = 1 (top) and n = 4 (bottom). 



able to reduce the error bar at least one order of magni- 
tude. Notice that for the Bertschinger-Zukin model, our 
error bars for an ideal experiment (no noise, Z max = 1000) 
compared with those in [41] (a = 0.61) are improved by 
a factor 3.74. This is due to the fact that in our analysis 
we account for the weak lensing in the power spectrum. 



B. Bispectrum 

As already mentioned in the previous section, even 
for perfectly Gaussian primordial perturbations ISW and 
WL effects generate some non-Gaussian features in the 
CMB. In this subsection, we explore the possibility to 
constrain cosmological parameters via the ISW-lensing 
bispectrum. For a rotation-invariant sky, all the cos- 
mological information is encoded in the angular-averaged 
bispectrum (e.g., see [58]), defined as: 



Bi. 



(2Ji + l)(2/ 2 + l)(2J 3 



47T 




i b hhh 

(30) 

where the reduced bispectrum bi r i 2 i 3 produced by the 
ISW-Lensing correlation is given by Eq. (27). The aver- 
aged bispctrum is different from only for li, I2, 13 triples 
that satisfy triangle inequality and parity conservation 
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i n=6, PLANCK 





{ 

1 


0.2 0.4 0.6 0.8 




1 


) 0.01 0.02 0.03 0.04 0.05 




0.8 


BZ model, ideal 


0.8 


BZ model, PLANCK 




likelihood 




likelihood 






u 

{ 


0.2 0.4 0.6 0.8 


u 


) 0.01 0.02 0.03 0.04 0.05 


FIG. 6. Hu-Sawicki model (n=6) and Bertschinger-Zukin 
model power spectrum likelihood. 


TABLE III. B error bars (68% CL) from bispectrum X 2 for 
the Hu-Sawicki model 


n 


WMAP 




PLANCK 




1 


5.42 x 10" 2 




2.90 x 10" 4 




4 


5.10 x 10" 1 




1.89 x 10" 3 




6 


5.64 x 10" 1 




6.29 x 10~ 3 





relationships: 



h < h < h , 

h + h + h = even , 

\lj — h\ < h < lj + h ■ 



(31a) 
(31b) 



Following the approach of the previous section, we now 
consider the chi-squared variable [37, 38, 58, 59]: 



Xb 



E 



B, 



fid 



where the bispectrum variance reads 

2 a r T r T r T i r T _ 

with the weights 



(32) 



Ni) , (33) 



(h =h = h) , 
(two identical) 
(all different) . 



(34a) 



The CMB bispectrum likelihood for the Hu-Sawicki 
models with n = 1,4,6 and the Bertschinger-Zukin 
model are plotted in Fig. 7 and Fig. 8, respectively, and 
the corresponding B error bars are listed in Tab. Ill and 
Tab. IV. Compared with the B error bars obtained from 
the power spectrum analysis, we do not find any signifi- 
cant improvement from the bispectrum alone. In details, 



TABLE IV. B error bars (68 ( 


% CL) from bispectrum x 2 f° r 


the Bertschinger-Zukin model 




Ideal (Z max = 1000) 


PLANCK 


1.92 x 10 _1 


3.85 x 10~ 2 



i n=1, PLANCK 



n=1,'wMAP 1 ^ ' — 

0.2 I 

— — ^^^^^^^^ 



0.02 0.04 0.06 0.0 



0.0002 0.0004 0.0006 0.0008 0.001 



2 0.4 





0.002 0.004 0.006 



FIG. 7. WMAP/PLANCK experiments bispectrum likeli- 
hood of Bq, for n = 1 (top) and n — 4 (bottom). 



for the Hu-Sawicki models with n — 4, 6 the ratio be- 
tween the error bars from the power spectrum and those 
from the bispectrum are 1.09 and 2.40, respectively. For 
the Hu-Sawicki model with n = 1 and the Bertschinger- 
Zukin model, the error bars from the bispectrum are even 
larger than those from the power spectrum. In the fol- 
lowing section we are however going to explore the possi- 
bility of improving our constraint through a joint analysis 
of power spectrum and bispectrum. 



C. Joint Analysis 

As already specified in the last part of Sec. Ill, the 
power spectrum can correlate with the ISW-Lensing 
bispectrum via both the primordial and the secondary 
anisotropies (ISW and weak lensing). However, compar- 
ing with the primordial and ISW channels, the cross- 
correlation through weak lensing is sub-dominant. The 
reason is that cross-correlator of the power spectrum and 
bispectrum via weak lensing channel is an 8-point func- 
tion of a Gaussian random fluctuations (a; m , (pi m ), while 
both the primordial and ISW are a 6-point function. 
Therefore, wc will focus on the unlcnscd cross-correlation 
channel in what follows, i.e. the power spectrum cross- 
correlation with the bispectrum through unlensed fluctu- 
ations. 

If wc compose the power spectrum and the bispectrum 
into a vector A = (AC/, Ai?; 1 ; 2 ; 3 ), the joint likelihood 
function reads 



£ oc cxp (A • CoV -1 • A T ) 



(35) 
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r 



i BZ model, PLANCK 



0.04 0.08 0.12 0.16 0.2 



FIG. 8. Hu-Sawicki model (n=6) and Bertschinger-Zukin 
model bispectrum likelihood. 



with 



c, 



E 



21 + 1 



^ T =E 



21 + 1 



(40) 
(41) 



We emphasize again that here we cross-correlate ISW- 
Lensing bispectrum with unlensed power spectrum. 
Prom the above result, we can see that the first two terms 
in Eq.(39), which peak at the squeezed configurations, 
give the main contributions to the cross-correlation, while 
the last term only affect the equilateral configurations. 
Since our ISW-Lensing bispectrum signal peaks at the 
squeezed limit, the contributions from the third term is 
negligible. In Appendix. B, we show a computation ex- 
ample about <7p B . 



where the inverse of the covariance matrix can be written 
as a 2 x 2 blocked form 



CoV" 1 



alal - (a 2 B ) 2 



2 

7 PB,iZiZ 2 i 3 



When the power spectrum and the bispectrum are weakly 
correlated (erp B <C fp B ), the above formula can be re- 
duced into [60] 



C cx exp(-x 2 /2) , X 2 = xl + xl + xIb ■ (37) 
For the computation of cross likelihood we have 



Xpb = 



E 



-2al BMihh x AC? x ASjfw-L 



where the variance reads 



ct 2 er 2 

°P,l u B,l 1 l 2 h 



(38) 



ee AC?AB}™- L ) = 



(2Ji + l)(2/ 2 + 1) (2/a + l) I h h h 

I 



4tt 

/l(/l + l)+/ 2 (/ 2 + D r / 3 (/ 3 + l) 

-L 



Cfct?CZ)+ Sperm.} - C[ B^ h 



(2/i + l)(2/ 2 + 1) (2/3 + 1) / h h h 

\ 



4tt 

-/l(/l + l)+/ 2 (/2 + l) + / 3 (/ 3 + l) 



(2/ 3 + 1) 

(2i + iy Cl Cl * Cl * 



(2l 2 + l) Cl C >- 



+ 5perm. 



(39) 




FIG. 9. Likelihood ratio between cross-correlation and bis- 
pectrum Xpb/Xb- 



In Fig. 9, we show the likelihood ratio between cross- 
correlation and bispectrum for several Bq values. This 
result indicates that the maximum correction from cross- 
correlations are 3%, so we can safely neglect the cross- 
correlations between the power spectrum and the bispec- 
trum in the joint analysis. In Fig. 10, we plot the joint 
likelihoods of Bo for the Planck experiment, and the cor- 
responding 68% CL results are listed in Tab.V. Compar- 
ing with the error bars from the power spectrum alone, 
there is an improvement of a factor ranging from 1.14 to 
5.32 (depending on the specific modified gravity model). 



V. CONCLUSIONS 

The correlation between the Integrated Sachs Wolfe 
(ISW) effect and gravitational lensing is expected to 
generate a bispectrum signature in the CMB . In fact 
the Planck satellite can detect this secondary non- 
Gaussianity [42-47] . The high precision measurements of 
CMB anisotropics expected from Planck data will enable 
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5e-05 0.0001 0.00015 0.0002 
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FIG. 10. Joint analysis results: spectrum and bispectrum 
product, no cross-correlation included. 



TABLE VI. B error bars (68% CL) summary for Planck 
experiment 



Model 


Pow 




Bisp 




Joint 




1 


6.68 x 10" 


5 


2.90 x 10" 


4 


5.85 x 10" 


-5 


4 


2.06 x 10" 


•3 


1.89 x 10" 


-3 


6.90 x 10 


-4 


6 


1.51 x 10" 


-2 


6.29 x 10" 


-3 


2.84 x 10" 


-3 


BZ 


1.10 x 10" 


-2 


3.85 x 10" 


-2 


7.87 x 10" 


-3 



small value of n and the Bertschinger-Zukin model. Fur- 
themore, we have provided the full formulae to compute 
the cross-correlation between the power spectrum and 
the bispectrum (that are necessary in the joint-analysis 
approach). One of our main results is that this cross- 
correlation can be safely neglected and this conclusion 
seems quite model independent. 



TABLE V. Bo error bars (68% CL) from the power spectrum- 
bispectrum joint analysis 



n/BZ Joint Power/Joint 

1 5.85 x 10~ 5 1.14 

4 6.90 x 10~ 4 2.99 

6 2.84 x 10~ 3 5.32 

BZ 7.87 x 10~ 3 1.40 
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us to detect such signature at a high level of significance 
[42-47]. This will open the interesting possibility of ac- 
tually using ISW-lcnsing bispectrum measurements as a 
tool to constrain cosmological parameters (see, e.g., [37- 
41]). Since the ISW-Lcnsing bispectrum is composed by 
the product of primordial power spectrum and the ISW- 
Lensing cross-spectrum, it carries information on both 
early (recombination) and late-time epochs. In this paper 
we investigated the possibility to use forthcoming Planck 
data to constrain some modified gravity models and their 
deviations from GR as a candidate to explain the late- 
time acceleration of the universe, focusing in particular 
on the possibility to perform a joint analysis of CMB 
power spectrum and ISW-Lensing bispectrum. 

Firstly we compared the ability to constrain cosmologi- 
cal parameters via CMB power spectrum only in WMAP 
and Planck missions. Due to the high resolution and sen- 
sitivity of the Planck dataset, Planck data will allow to 
improve the constraint on the modified gravity param- 
eter B at least by one order magnitude compared to 
what one can do with WMAP. Secondly we performed a 
joint power spectrum/bispectrum analysis. Our results 
show that the joint analysis could reduce the B 1-a er- 
ror bar by a factor ranging from 1.14 to 5.32 w.r.t. to 
what can be achieved with the power spectrum alone. 
The improvement factor depends on the model consid- 
ered: for the Hu-Sawicki models with large value of n 
the improvements are quite significant unlike those with 



Appendix A: CMB power spectrum Fisher matrix 
analysis 

In this Appendix, we present the results of our Ci 
Fisher matrix analysis. The Fisher matrix for CMB 
anisotropics and polarization is [50-55] 



- 1 da 



Yl 



I X.Y 

where the covariance matrix is defined as 



XY dpj 



Cov 



tt {21 + l)/ sky 



(Cov ; ) 
(Cov ; ) = 

V / TE 



EE (21 + l)/ sky 

2 



Cti + w^b; 



Cei+w^B- 



TE (21 + l)/ sky 



r 2 

— ' f 



CI ! 



with 



, v -2 

W X ,c = {vX.cObcam) , % = )^ W X ,, 

c 

-V»xAi, Xe(T,P) 



B, 



B, 



w x 
~ exp 



8 In 2 



(Al) 

(A2) 
(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
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TABLE VII. WMAP and Planck experiment 



Experiment 


Frequency 


^bcam 


<7t 


Op 


WMAP: 


94 


12.6 


49.9 


70.7 




60 


21.0 


30.0 


42.6 




40 


28.2 


17.2 


24.4 


Planck: 


217 


5.0 


13.1 


26.8 




143 


7.1 


6.0 


11.5 


Junk:2012qt 


100 


10.0 


6.8 


10.9 




70 


14.0 


12.8 


18.3 



arcminutes. Sensitivities <jt and op are in fiK per FWHM 
beam. 



In the above expressions, Bf is the Gaussian beam win- 
dow function ,#bcam is the beam's full-width at half- 
maximum (FWHM), and ax are the sensitivities. In this 
paper we choose beams and sensitivities representative of 
WMAP and Planck experiments, following specifications 
from [55, 57]. For convenience we list channels, beams 
and noise levels we used in Tab. VII. 

As shown in the main text, the likelihood of Bq is 
highly non-Gaussian, which makes forecasting B error 
bars through a Fisher matrix analysis unfeasible. On the 
other hand we still use the Fisher matrix as a way to get a 
reasonable order of magnitude estimate of the correlation 
between different parameters. 

We jointly analyze 7 cosmological parameters within 
the Hu-Sawicki model (A s ,n s ,T, f2b/i 2 , Q m h 2 , h, B ). The 
first 6 are the standard ACDM parameters, while Bo is 
the new modified gravity parameter. Actually, as stated 
in Sec. II B the Hu-Sawicki model is charaterized by 2 
parameters (Bo,n), however, in our calculations we fix 
n for a few specific values. The are various reasons for 
this choice: firstly, there exists a parameter degeneracy 
between n and B ; secondly, the current constraints on 
the parameter n is rather weak [9] and the results are 
strongly dependent on the value of n. Given the above 
considerations, we fix n for a few specific values in the 
following numerical analysis, i.e. we can deal with the 
Hu-Sawicki model with different n's as different models. 

Our Fisher matrix results arc summarized in Ta- 
ble VIII. The first six rows are obtained for the Hu- 
Sawicki model with n = 4. Actually, we find that the 
dependence of these standard parameters on n is rather 
weak, i.e. we can obtain similar results for other values of 
ij. As a numerical check, we verified that our results for 
the 6 "standard" parameters match well both WMAP- 
7yr [56] and Planck blue book results [57]. We find no 
significant degeneracies between B and other parame- 
ter, thus justifying the simple 1-parameter likelihood ap- 
proach we use in the main text to forecast Bq error bars 

Furthermore, in Tab. IX we list the cross-correlation 
coefficients of several cosmological parameters with B in 
Hu-Sawicki model with n = 4. Our results show that all 
the cross-correlation coefficients are small. 



TABLE VIII. Fisher matrix joint analysis results (68%CL) on 
lensed power spectrum for the Hu-Sawicki model 



/ ' I > 


Fiducial value 3. 


WMAr 




r .LAIN Urv 


A 

A s 


2.43 X 1U 


nf?i w in — 

y.bi X 1U 


ii 


o ac w in — 11 
3.4b X 1U 


n s 


0.963 


1.15 x 10" 


-2 


3.38 x 10~ 3 


T 


0.088 


1.58 x 10" 


-2 


3.92 x 10~ 3 


Q b h 2 


0.0226 


5.53 x 10" 


-4 


1.47 x 10~ 4 




0.1109 


5.31 x 10" 


-3 


1.23 x 10 -3 


h 


0.70 


1.72 x 10" 


-2 


6.20 x 10~ 3 


B (n = 4) b 





1.81 x 10" 


-3 


8.33 x 10~ 5 


B (n = 1) 





4.41 x 10" 


-4 


1.92 x 10~ 5 



a Fiducial values of the first six parameters are taken as 

WMAP7yr best-fit values. 
b Results of the first six parameters include polarization 

contributions, while those of Bq are only from Cf T mode. 



TABLE IX. Cross-correlation coefficients of the standard cos- 
mological parameters with Bo in Hu-Sawicki model (n = 4) 



Bo-CP 



cross- correlation coefficients 



n s 

T 

Q b h 2 

h 



-7.9 x 10~ 2 

-0.12 
-8.0 x 10~ 2 

-0.18 
-4.7 x 10~ 2 
8.2 x 10~ 2 



Appendix B: Computation example of op B 

As an example of the computation for the cross- 
likelihood, Eq. 39 we show explicitly the following term: 



Cfcf^Cl) - CfBf 



jISW-L 



E 



(2/ + l)(2/ 2 + l)(2/ 3 + l) 

•T dISW-L 

(2/ + l)(2Z 2 + l)(2/ 3 + l) 

* T 



m ' , m 2 , m 3 , rn , m 2 , m 3 
a lm a lm<Pl2m2 a h rn 3 I (1(2(3 



E 

m / ,m 2 ,m3,m,m2,m3 

1 



1 



,T* 



+ 2 \ a itn a i pi^ 2 m 2 a hf ri2 a hm 3 a hm 3 



F* T 1* T T* T 

: |m Ct ip ( ? ) i2 | m,2 a l2jTl2 a l3^l3 a h^ l3 



(Bl) 
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where we have subtracted 3 disconnected terms. 
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